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Abstract: 

Zeilberger's algorithm provides a method to compute recurrence and differential equations from 
given hypergeometric series representations, and an adaption of Almquist and Zeilberger com- 
putes recurrence and differential equations for hyperexponential integrals. Further versions of 
this algorithm allow the computation of recurrence and differential equations from Rodrigues 
type formulas and from generating functions. In particular, these algorithms can be used to 
compute the differential/ difference and recurrence equations for the classical continuous and 
discrete orthogonal polynomials from their hypergeometric representations, and from their Ro- 
drigues representations and generating functions. 

In recent work, we used an explicit formula for the recurrence equation of families of classical 
continuous and discrete orthogonal polynomials, in terms of the coefficients of their differen- 
tial/difference equations, to give an algorithm to identify the polynomial system from a given 
recurrence equation. 

In this article we extend these results be presenting a collection of algorithms with which any of 
the conversions between the differential/ difference equation, the hypergeometric representation, 
and the recurrence equation is possible. 

The main technique is again to use explicit formulas for structural identities of the given poly- 
nomial systems. 

1 Structural Formulas for Classical Orthogonal Polynomials 

A family 

y{x)=p n {x) = k n x n + ... (n£N := {0,1,2,... },k n / 0) (1) 

of polynomials of degree exactly n is a family of classical continuous orthogonal polynomials 
if it is the solution of a differential equation of the type 

a{x) y"{x) + t{x) y'{x) + \ n y{x) = (2) 

where cr(x) = ax 2 + bx + c is a polynomial of at most second order and t(x) = dx + e is a 
polynomial of first order. Since one demands that p n (x) has exact degree n, by equating the 
highest coefficients of x n in (2) one gets 

X n = —(an(n — 1) + dn) . (3) 

Similarly a family p n (x) of polynomials of degree exactly n, given by (1), is a family of discrete 
classical orthogonal polynomials if it is the solution of a difference equation of the type 

a{x) AVy(x) + t{x) Ay(x) + \ n y{x) = , (4) 
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where 

Ay{x) = y(x + 1) - y(x) and Vy{x) = y{x) - y(x - 1) 

denote the forward and backward difference operators, respectively, and o{x) = ax 2 + bx + c 
and t{x) = dx + e are again polynomials of at most second and of first order, respectively. 
Again, (3) follows. 

Since AV = A — V, (4) can also be written in the equivalent form 

(a{x) + r{x))Ay{x) - a{x) Vy{x) + \ n y{x) = , 

and replacing x by x + 1 we arrive at 

(g{x + 1) + t{x + 1)) A 2 y{x) + t(x + 1) Ay(x) + X n y(x + 1) = . (5) 

It can be shown (see e.g. [14] or [16]) that any solution p n (x) of either (2) or (4) satisfies a 
recurrence equation 

p n+ i(x) = (A n x + B n )p n (x) -C n p n -i(x) (n€N ,p_i = 0) (6) 

or equivalent ly 

xpn(x) = a n p n+ i(x) + b n p n (x) + c n p n _i(x) (7) 

with 

1 , B n C n . . 

an = -7- , K = —7- C n = — . (8) 

-n-n ^-n 

In [14] (compare [16], Eqs. (5) and (10)) we showed that the coefficients A ni B ni and C n are 
given by the explicit formulas 

_ K+i 

2b n {an + d — a) — e (—d + 2 a) k n+ \ 
(d + 2an)(d — 2a + 2an) k n 

and 

{{an+d— 2a)n(Aca — b 2 )+4a 2 c— ab 2 +ae 2 —4acd+db 2 — bed+d 2 cj {an+d—2a)n fc n+l 
(d—2a+2an) 2 {2an — ?>a+d) {2an — a+d) k n -\ 

in the continuous case, and by the formulas 

_ k n+ i 

n(d + 2b) (d + an — a) + e(d — 2a) k n+ i 
( 2 an — 2 a + d) ( d + 2 an ) k n 

and 

C n = — ^( n — 1 ) ( d+an — a) ( a n d—db—ad+a 2 n 2 — 2a 2 n+4ca+a 2 + 2ea — b 2 ) — dbe+d 2 c+a e 2 ^ 

{an + d— 2a) n k n+ \ 
{d—a + 2an) {d+2an — 3a) {2an — 2a + d) 2 k n -\ 
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in the discrete case, in terms of the coefficients a, 6, c, d and e of the given differential/ difference 
equation. 

Orthogonal polynomials satisfy further structure equations. One of those is given by the 
derivative/ difference rules (see e.g. [14]) 

a{x)p' n {x) = a n p n+ i(x) + P n Pn(x) +-y n Pn-i{x) (n G N := {1,2,3,...}) , (9) 

and 

<*{x) Vp n {x) = a n p n+ i{x) + PnPn(x) + ^ n Pn-l{x) (n € N) , (10) 



or 



^cr(x) + t(x)^ Ap n (x) = S n p n+ i(x) + T n p n (x) + i? n p n _i(x) (n G N) , (11) 
respectively. Here 

S n = 0> n , T n = (3 n — X n , -Rn = 7n • (12) 

In [14] we showed that the coefficients a n ,/3 n , and 7 n are given by the explicit formulas 



= an- 



n (a n + d — a) (2 e a — db) 



(d + 2an)(d-2a + 2an) 
and 

((n — l)(an+d— a)(4ca — 6 2 )+ae 2 +d 2 c— &ed)(an+d— a)(an+d— 2a)n k n 



In 



(d—2a + 2an) 2 (2an — 3a + d) (2an — a + d) k n -\ 
in the continuous case, and by the formulas 



an = an 



n ( d + an — a) (2 an d — a d — db + 2 e a — 2 a 2 n + 2 a? n 2 ) 
^™ ~~ (2an - 2a + d) (d + 2an) 

and 

7 n = ((n — l)(d+an — a)(and— db — ad+a 2 n 2 — 2a 2 n+Aca + a 2 +2ea — b 2 ) — dbe + d 2 c+ae 2 ^j 

(d+an — a)(an + d—2a)n k n 
(d—a + 2an) (d+2an — 3a) (2an — 2a + d) 2 k n -\ 

in the discrete case, respectively. 

Now, we develop further structural identities. Taking the derivative in (2), we get 

= o{x)p"{x) + (t{x) + a'(xfj p'n{x) + (a„ + t'{x)) p' n (x) 

= (ax 2 + bx + c)p'n{x) + ((d + 2a) x + (e + &)) p'^x) + (A„ + d)p' n {x) , 

hence y(:r) :=p^(:r) satisfies a differential equation 

(a'x 2 + b'x + c') y"(a;) + (d'x + e') y'(a;) + \' n y{x) 
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of the same type with 

a' = a , b' = b , c' = c , d' = d + 2a , and e' = e + 6 . (13) 
Prom this we deduce that the equation 

xp'Jx) = a^p' n+1 {x)+P*p' n {x)+'y*p' n _ 1 {x) , (14) 
namely a recurrence equation for p' n (x), is valid, and from (13) it follows that 
a* = a n (a, b,c,d + 2a, e + b) , /?* = 6 n (a, 6, c, d + 2a, e + b) , 

and 

7* = c n (a, b,c,d + 2a, e + b) , 

where a n (a,b,c,d,e),b n (a,b,c,d,e), and c„(a, 6, c, d, e), are given by (8) and the explicit for- 
mulas for A n ,B n and C n . 

Similarly in the discrete case, applying A to (4), we get for y(x) := Ap n (x) 

= (a(x + 1) - Aa(x) ) A\7ij(x) + (r(a; + 1) + Aa(xf) Ay(x) + (A„ + Ar(s)) y(a;) 
= (aa; 2 + 6a; + c) AVy(a:) + ((d + 2a) a; + d + e + a + ft) Ay(x) + (A„ + cZ) y(a;) , 

hence y(a;) := Ap n (x) satisfies a difference equation 

(ax 2 + b'x + c ) AVy(x) + (d'x + e ) Ay(x) + A'„ y(a;) 
of the same type with 

a' = a , 6' = b , c' = c , d' = d + 2a, and e' = d + e + a + 6. (15) 
From this we deduce that the equation 

a: Ap„(a:) = a* Ap„+i(a:) + /3* Ap n (x) + 7* Ap„_i(a:) , (16) 
namely a recurrence equation for Ap n (x), is valid, and from (15) it follows that 

a* = a n (a, b,c,d + 2a, d + e + a + b) , /3* = 6 n (a, 6, c, d + 2a, (Z + e + a + 6) , 

and 

7n = c„ (a, 6, c, d + 2a, d + e + a + 6) , 

where a n (a,b,c,d,e),b n (a,b,c,d,e), and c n (a,b,c,d,e), are given by (8) and the explicit for- 
mulas for A n ,B n and C n . 

To obtain a derivative rule for y(x) := p' n {x), we take the derivative of (9) to get 

a{x)p„{x) +a'(x)p' n (x) = a n p' n+1 (x) + (3 n p' n (x) +-y n p' n _ 1 (x) . 
Applying (14) to replace xp' n (x) results in a derivative rule of the form 

a{x)p„{x) = d n p' n+l {x) + b' n p' n {x) + c n p' n _ x (x) . (17) 
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Similarly in the discrete case a difference rule of the form 

a{x) AVp n (x) = a' n Ap n+l (x) + b' n Ap n (x) + c' n Ap^x) (18) 
can be obtained for y(x) := Ay(x). 

Finally we substitute (17) in the differential equation. This gives 

a 'nPn+l( x ) + b 'nPn( x ) + c 'nPn-l( x ) + T ( x ) Pn( x ) + K.Pn(x) = , 

and replacing xp' n (x) by (14), again, we obtain an equation of the type 

Pn{x) = a n p' n+1 (x) + b n p' n {x) + c n p' n _ 1 (x) , (19) 

in the continuous case, and a similar procedure gives 

p n {x) = a n Ap n+ i(x) + b n Ap n (x) + c n Ap n _i(x) (20) 

in the discrete case. Note that in the discrete case also corresponding equations concerning V 
are valid. 

We note in passing that our development shows by simple algebraic arguments that whenever 
p n { x ) is a polynomial system of degree exactly n, satisfying a differential/difference equation 
of type (2)/(4), a recurrence equation of type (6) and a derivative/difference rule of type 
(9) / (10) then the system p' n+1 (x) (Ap n+ \(x)) is again such a system. This has nothing to do 
with orthogonality. Indeed, in our further development it will become rather important that 
in the continuous case the powers x n and in the discrete case the falling factorials 

x R := x(x -!)■■■ (x-n + 1) = (x-n + l) n = (-1)" (-x) n 

which by no means become orthogonal families, have these properties. 

To deduce the coefficients a* , (3*, 7*, a' n , b' n , c' n , and a n ,b n ,c n , we can follow the above instruc- 
tions, or we apply the following method: Substituting 

Pn(%^) ~ kri X ~\~ k n X -\- k n X 

in the differential/ difference equation and equating the coefficients of x n determines A n , while 
equating the coefficients of x n ~ 1 and x n ~ 2 gives k' n , and respectively, in terms of k n . These 
values can be substituted in p n (x). Next, we substitute p n (x) in the proposed equation, and 
equate again the three highest coefficients successively to get the three unknowns in terms of 
a, 6, c, d, e, n, k n -\ , k n , and k n+ \ by linear algebra. 

These computations can be easily carried out by a computer algebra system, e.g. by Maple. 
With few seconds of computation time, we get 

Theorem 1 For the solutions of (2) and (4), the relations (14), (17), (19), and (16), (18), 
(20), respectively, are valid. The coefficients «*,/?*, 7*, a' n ,b' n ,c' n , and a n ,b n ,c n , are given by 

* _ n kn 

Un ~ n + l ' k^l ' 

* — 2bn(an + d — a)+d{b — e) 
^ n ~ (d + 2an){d-2a + 2an) ' 
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„ ((n — 1) (an + d — a) (4ca — b 2 ) + ae 2 + d 2 c — bed) n(an + d — a) 

ln = 1 " ' 



(d — 2a + 2an) 2 (2an — 3a + d) (2 an — a + d) k n -\ 
, an(n — 1) k n 



n + l k 



n+l 



, (n — 1) (an + d) (2 e a — db) 
n ~~ (d + 2an)(d-2a + 2an) 



i 

c n 



( ( n — 1) (an + d — a) (4c a — b 2 ) + ae 2 + d 2 c — bed) (an + d) (an + d — a)n k r , 



(d — 2a + 2 an) 2 (2 an — 3 a + d) (2 an — a + d) k n -\ 
^ 1 k n 



n + l k n+ i 
2ea — db 



(d + 2an)(d-2a + 2an) 
( ( n — 1) (an + d — a) (Ac a — b 2 ) + ae 2 + d 2 c — bed) an k r , 



( d — 2 a + 2 a n ) 2 ( 2 an — 3 a + d) ( 2 an — a + d) k n -\ 
in the continuous case, and 



0* 
r-'n 



n + 1 k n+1 
-n(d + 2a + 2b) (d + an — a) — d(e — a — b) 
(2an - 2a + d) (d + 2an) 



7* = — ((n — 1) (d+an — a) (and—db—ad+a 2 n 2 —2a 2 n+4:ca + a 2 +2ea — b 2 )—dbe + d 2 c+ae 2> j 

(d + an — a)n k n 
(d — a + 2an) (d + 2an — 3a) (2an — 2a + d) 2 k n -\ ' 

, an(n — 1) k n 
Gn ~ n + l ' k^[ ' 

(n — 1) (an + d) (2and — ad — db + 2ea — 2 a 2 n + 2a 2 n 2 ) 



b' 



(2an-2a + d)(d + 2an) 



c' n = (^(n — l)(d+an — a)(and—db—ad+a 2 n 2 — 2a 2 n+4:ca+a 2 + 2ea — b 2 ) — dbe+d 2 c+ae 2 ^ ■ 

(d + an — a) (an + d) n k n 
(d — a + 2 an) (d + 2an — 3 a) (2an — 2 a + d) 2 k n -\ ' 

a — ^ ^ n 
n + l k n+ i 

^ — 2an(d + an — a) — db + ad — d 2 + 2ea 
n ~ (2an-2a + d)(d + 2an) ' 

c n = {(n—1) (d+an — a) (and— db — ad+a 2 n 2 — 2a 2 n + 4ca + a 2 + 2ea — b 2 )—dbe + d 2 c+ae 2 ^j ■ 

an k n 
(d — a + 2an) (d + 2an — 3a) (2an — 2a + d) 2 k n -\ 

in the discrete case. □ 
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Note that (19) gives an immediate formula for the antiderivative of a continuous orthogonal 
polynomial in terms of its neighbors, so that definite integrals can easily be computed, whereas 
(20) gives an immediate formula for the antidifference of a discrete orthogonal polynomial in 
terms of its neighbors, so that definite sums can easily be computed. 

As a direct consequence of Theorem 1 we have the following representations. The definition 
of the continuous and discrete families will be given in § 3 and § 5. 

Corollary 1 The classical continuous orthogonal polynomials have the following antideriva- 
tive representations: 

J H n (x) dx = - — - - H n+1 (x) 

(see e.g. [24], (5.5.10)), 

j L^(x)dx = L^(x)-LH l (x) 

(see e.g. [25], VI (1.14)), 

[BPWdx = 7 lWo 2(n + 1 + " } 

J n v ' (n + 1) (2n + a + 1) (2n + a + 2) n+n ' 

4 , n 

-i BS a >(x\ 

(2n + a)(2n + a + 2) n { ' 

, R («) / \ 

(n + a)(2n + a)(2n + a+l) n ~ l[X) ' 

f 1 1 

/ C${x) dx = — — — C% +1 (x) dx - — — — C"_i(a;) dx 
J 2 [n + a) 2 [n + a) 

(see e.g. [25], V (7.15)), 

J n () ~ {2n + a + f5+l){2n + a + f5 + 2) Fn+l[X)dX 

+ P («,/3) (*) ^ 

(2n + a + (3) (2n + « + /? + 2) n v y 

2 ( n + g) ( n + g) P («:/3) M rf _ 

(n + « + /3)(2n + a + /3)(2n + a + /3+l) n " x lJ 

(see [12], Theorem 6). 

The classical discrete orthogonal polynomials have the following antidifference representations: 

(x) = — c ^n+l( x ) , 

x 

5>^(z,7V) = A^fo JV) -pfc^^JV) , 



X 



1 1 (n — N) (n + TV) 

^ iV ) = 2(2^TT) iV) - 2 ^ 7V) + 2(2n + l) *° 
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n + a + f3 + 1 



(2n + a + (3 + 1) (2n + a + /3 + 2) n+1 1 ' ; 

2n 2 + 2n + 2na + 2n/3 + a - aN + /37V + a(3 + /3 + /3 2 



(2n + a + f3) (2n + a + /3 + 2) 
, (n + a) (n + g) (n - N) (n + a + g + 7V) (aJ) 

(n + a + /3)(2n + a + /3)(2n + a + /3+l) "- 1 1 ' >' 

Proof: Using the representations for a n , b n and c n of Theorem 1 with the particular values 
for a,b,c,d,e and k n of the families (see e.g. [1], [19]) give the results. □ 

Note that the representations for a n , b n and c n of Theorem 1, if applied to P n (x) = x n or 
Pn{x) = x—, respectively, yields the simple results 



/ 



x n dx 



n+1 



„«+i 



and 



V iE a = _!_ a ^±l j 
respectively. The latter is equivalent to the well-known identity 



£(T) 



k=0 



1 m i m+n -• 



A;=0 



A;=n 

(n + 1)! V m ; 



k=m+n 
k—n 



The polynomial system 



-n 



-x — n- 



1+/3 



1 



: (-l) n - 2j Po 



-n, x + 



a , (21) 



which was given in [14], is not orthogonal, but Theorem 1 is still applicable, and we get 

a (n + 1) n+1 



2 Connection Coefficients 



In this section we would like to consider the problem to determine connection coefficients 
between different polynomial systems. Here we assume that P n (x) = k n x n + ... (n 6 No) 
denotes a family of polynomials of degree exactly n and Q m (x) = k m x m + . . . (n 6 No) denotes 
a family of polynomials of degree exactly m. Then for any n £ No a relation of the type 

n 

P n (x)=J2C m (n)Q m (x), (22) 

is valid, and the coefficients C m (n) (n £ No, m = 0, ...,n) are called the connection coef- 
ficients between the systems P n (x) and Q m (x). For simplicity we assume that C m (n) are 
defined for all integers n,m and that C m (n) = outside the above n x m-region. 
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The connection coefficients between many of the classical orthogonal polynomial systems had 
been determined by different kinds of methods (see e.g. [24], [10], [20]) until Askey and Gasper 
[6] used recurrence equations to prove the positivity of the connection coefficients between 
certain instances of the Jacobi polynomials. In a series of papers ([21]-[22], [3]) Ronveaux 
et. al. recently used such a method more systematically. Here we will present an algorithmic 
approach different from theirs. 

Hence, the main idea is to determine recurrence equations for C m (n). Since C m (n) depends on 
two parameters m and n, many mixed recurrence equations are valid as we shall see. The most 
interesting recurrence equations are those which leave one of the parameters fixed. We will 
determine those recurrence equations, hence pure recurrence equations w.r.t. m and n. The 
success of this method will heavily depend on whether or not these recurrence equations are of 
lowest order, i.e., whether or not no recurrence equations of lower order for C m (n) are valid. In 
cases when the order of the resulting recurrence equation is one, it defines a hypergeometric 
term which can be given explicitly in terms of shifted factorials (or Pochhammer symbols) 
(a)k = a(a + 1) ■ ■ ■ (a + k — 1) = T(a + k)/T(a) using the initial value C n (n) = k n /k n . We will 
see that there are many instances for this situation. 

Note that w.l.o.g. we could assume that k n = k m = 1, i.e., that P n (x) and Q m (x) are monic 
systems with connection coefficients C m (n), because if P n (x) and Q m (x) have leading coeffi- 
cients k n and k m , respectively, then their connection coefficients C m (n) are given by 

C m (n) = y~ C m (n) . 

In the last section we have already solved a rather special connection problem: (19)/(20) ex- 
presses the connection between the polynomial systems P n (x) = p n (x) and Q m (x) = p' m+ i(x) 
or Q m (x) = Ap m+ i(x), respectively. In this case the connection coefficients turn out to be 
rather simple: almost all of them (namely all with m < n — 2) are zero. 

Now, we consider the generic case. We assume that P n (x) is a polynomial system given by 
(2)/(4) with a(x) = ax 2 + bx + c, and t(x) = dx + e, and that Q m {x) is a polynomial system 
given by (2) / (4) with a(x) = ax 2 + bx + c, and f(x) = dx + e. 

We know then that both P n (x) and Q m (x) satisfy a recurrence equation (7) whose coefficients 
a n (a, 6, c, d, e), b n (a, 6, c, d, e), and c„(a, 6, c, d, e) were given explicitly in the last section. Note 
that we will denote all coefficients connected with Q m (x) by dashes. Hence we have 

P n+1 (x)+b n P n (x) + 

and 

Qm (^0 = "m Qm+1 (^0 Qm(^) Cm Qm—1 {%) j 

all of a n ,b n , c n ,a m ,b m ,c m given explicitly. 

In three steps, we will now derive three independent recurrence equations for C m (n). First we 
consider the term x P n (x) (see e.g. [22]). Using the defining equation of C m (n), and the two 
recurrence equations for P n (x) and Q m (x), we get 

xP n (x) = a n P n+ i(x) + (x) + c n P n -i(x) 

n 

= ^2 ( a n C m (n + 1) Q m {x) + b n C m (n) Q m {x) + c n C m (n - 1) Q m {x)^ 

n 

= ^2 C m (n)xQ m (x) 

m=0 
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Q m+ l{x) + b m Q m (x) + C m Q m -l(x)j . 

By appropriate index shifts, we can equate the coefficient of Q m {x) to get the "cross rule" 

a n C m (n + l) + b n C m (n) + c n C m (n- 1) = a m _i C m _i(n) + 6 m C m (n) +c m+ i C m+ i(n) . (23) 

To deduce a second cross rule in terms of the same variables C m (n + 1), C TO (n), C m (n — 1), 
C m _i(n) and C m +i(n), we examine the term x P^(x) (or iAP n (i) in the discrete case). 
Using both recurrence equations for the derivatives/differences 

x P' n (x) = a* n P' n+ iix) + (3* n P' n (x) + Y n PLi(x) 

and 

x Q m (x) = a* m Q' m+1 (x) + p m Q' m (x) + f m Q' m _ x (x) 

(or analogously 

x AP n (x) = a* n AP n+l (x) + 0* AP n (x) + 7* AP B _i(i) 

and 

xAQ m (x) = a* m AQ m+1 (x) + (3* m AQ m (x) + ^* m AQ m ^ 1 (x) 
in the discrete case), we get 

x P n (x) = a n P n+ \ (x) + f3 n P n (x) + 7 n P n _i (x) 

n 

= (<C m (n + l)Q' m (x)+P*C m (n)Q' m (x) +1 *C m (n-l)Q' m (x)) 

m=0 

n 

= C "i(n) x Q' m (x) 

m=0 

n 

= c m(n) (a* m Q' m+1 (x) +]3* m Q' m (x) + "Jm Q' m -i( x )) ■ 

Again, by appropriate index shifts, we can equate the coefficient of Q m {x) to get the cross 
rule 

a* n C m (n + 1) C m {n) +j* C m {n- 1) = a*^ C m _i(n) +~0* m C m {n) +7™+i C m+1 {n) (24) 
(and the same result in the discrete case). In a similar way the cross rule 

a n C m (n + 1) + b n C m (n) + c„ C m (n - 1) = a m _i C m _i(n) + b m C m (n) + c m+ i C m+ i(n) (25) 

can be obtained. It turns out, however, that this relation is linearly dependent from (23) and 
(24), and hence does not yield new information. 

Now, we specialize a little. First, we consider the continuous case. To obtain reasonably 
simple results, we assume furthermore that a(x) = a(x). We consider the term a(x) P^(x). 
Then, using both derivative rules 

a(x) P' n {x) = a n P n+ i(x) + (3 n P n (x) + 7„ P n -\{x) 
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and 

a(x) Q' m ( x) — OL m 

Qm+1 (x) + P m Qm (x) + 7 m Q m -l (x) , 

we get 

a(x)P' n (x) = a n P n+ i(x) + (3 n P n (x) + -f n P n - 1 (x) 
n 

= ^2 (a n C m (n + 1) Q m (x) + (3 n C m (n) Q m (x) + j n C m (n - 1) Q m {x)) 

n 

= Yj C m {n)o(x)Q' m (x) 

n 

Q m +l(x) + (3 m Q m (x) + 7 m Q m -l{xfj . 

Again, by appropriate index shifts, this results in the cross rule 

a n C m (n + l) + P n C m (n)+'j n C m (n-l) = a TO _i C m _i(n) + /3 m C TO (n) + 7 m+1 C m+1 (n) . (26) 

To obtain a pure recurrence equation w.r.t. m, from the three cross rules (23), (24), and (26) 
by linear algebra we eliminate the variables C m (n + 1) and C m (n — 1), and to obtain a pure 
recurrence equation w.r.t. n, we eliminate the variables C m _i(n) and C m +i(n). For simplicity 
we consider the monic case. 

Theorem 2 Let P n (x) be a monic polynomial system given by the differential equation (2) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e, and let Q m (x) be a monic polynomial system 
given by (2) with cf(x) = a(x), and t(x) = dx + e. Then the relation (22) is valid, C m (n) 
satisfying the recurrence equation 

— ( m — n) (am + d — a + an) (d + 2 am) ( d + a + 2 a m ) (d + 3 a + 2 am) 

( d + 2 a m + 2a) 2 C m (n) + (-dbnd + 2d am b + dbd + 2damb + 2dena 
+ dde + 2ddbm — mbd 2 — ed 2 —4a 2 m 2 e — m 2 abd + bnda — 2eda 
— 4a me — 4edam + 2m a e + 2ea n —2ea n — mabd + 2mdea 
+ 2me a 2 — bn 2 da)(d + 2am + 2 a) (rn + 1) (d + a + 2am) (d + 3 a + 2am) 
C m+ i (n) — (d + 2 am) (m + 1) ( —a m — 2 a + a n — d + d) ( a m + a n + a + d) 
(a b 2 m 2 — 4 a 2 m 2 c — 8 a 2 mc + 2amb 2 — 4admc + mb 2 d — 4a d c — ae 2 + ab 2 — cd 2 
+ bed - 4a 2 c + b 2 d)(m + 2) C m+2 {n) = 

with respect to m, with initial values C n (n) = l,C„+i(n) = 0. Furthermore the recurrence 
equation 

-(d + 2an) 2 (d - a + 2an) (d + 2an + 2a) (d + a + 2an) (-m + n + 2) 

(d + am + a + an) C rn (n + 2) + (d — a + 2an) (d + a + 2an) (n + 2) (d + 2an) 
(— 2ead — bd 2 n — 2ma 2 e + 2 m 2 a 2 e + 2a 2 en — 4e a 2 n 2 + 2 db d — 2dea 

O - O -000 

+ mbda — bd + 2eda — 4eadn + 2edan — abdn —bdna + 2a en —m abd 

— 4e a 2 n — dmb d + 2danb + 2 dan 2 b + 2ddbn + 2dmea 

— ed 2 + dde)C m (n + \) + (d + 2an + 2a) (n + 2) (n + \ ) (an — am + d — d) 

(an + am — a + d) (bed — ae 2 — d 2 c — 4 a end — 4 a 2 cn 2 + ab 2 n 2 + nb 2 d) C m (n) = 

with respect to n is valid. 
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Proof: Using the explicit representations given in the last section in combination with (23), 
(24), and (26), and elimination of C m (n+1) and C m (n— 1), or C m _i(n) and C m+ i(n), respec- 
tively, yields the results. □ 

Note that the recurrence equation given in Theorem 2 reduces to two terms, and hence can 
be represented by hypergeometric terms, for the connection between Laguerre polynomials 
(P n (x) = L^\x),Q m (x) = Lm\x)), and between the Gegenbauer polynomials (P n (x) = 
C%(x),Q m (x) = C^(x)). We will consider these and more cases by another method in § 4. 

Now, let's switch to the discrete case. There are two possibilities to obtain a relation similar 
to (26). Replacing the derivative by V, the same argument gives (26), again, valid for a(x) = 
a(x). If <j{x) + t{x) = cf(x) + t(x), we can replace the derivative by A, and adopt the above 
argument to get the relation 

S n C m (n+l)+T n C m (n) + R n C m (n-l) = S m ^ 1 C m ^ 1 (n)+T m C m (n) + R m+1 C m+1 (n) . (27) 

Hence we get 

Theorem 3 Let P n (x) be a monic polynomial system given by the difference equation (4) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e, and let Q m (x) be a monic polynomial system 
given by (4) with o(x) = <j(x), and t(x) = dx + e. Then the relation (22) is valid, C m (n) 
satisfying the recurrence equation 

( d + 2 am + 2a) 2 (d + 3a + 2 am) (d + a + 2 am) (d + 2 am) ( —m + n ) 

(an — a + d + am) C m (n) — (d + 2 am + 2 a) (d + 3 a + 2am) (d + a + 2am) (m + 1) 

— o — 2 — o o — 2 q o q o — 

(—add + an d + 2ead — 2a em +bd m — 2na m —2a nm — a nd 

o o — — 2 — o — o - q o o — 

+ a n d — and — 2eadm + 4ea m + badm +2a m n + adbm 

o — o ^? ^? — 2 o o — 'x — 2 o o o — — 2 

+ 2a mdn —4a m —2ad m —4a dm +dnd —2a m d — a md — d ma 

— 2a 2 em + 2a 3 mn 2 +andd — 2a A m 2 — 2a 2 en 2 — d 2 d— dbd — 2 a 3 m 4 — 2 a 2 m d 

— oo — o o — 

+ 2dmand + 2a m nd + 4 a em + 4 a em d — 2 ae n d + 2 a ndm — 2a mdn 

— — o — o — o 

— 3amdd — dm ad — 2dbmd — 2am bd — 2ambd — 5dm a +2ea n 

— ded — d 2 md + ed 2 + dbnd + an 2 bd — anbd) C m+ \(n) + ( m + 1 ) 

(d + 2am) ((i + am + a + on) (m + 2)(4a 2 c m 2 + 2 a 2 e m 2 — bd 2 m — b 2 a m 2 
+ 2e a dm — bad m 2 + 4dcam + 2 a 2 c? + 4 a 3 m — 2 adbm + 4 a 3 m 3 

— o — 2 o o — *3 o — — 2 o ^ o 

— d6 m + ad m +2a dm +6a md + 2d ma + 4a em + 6a m 

+ 4dca-2b 2 am + 8a 2 cm + 2a 2 e-b 2 a + a 3 m 4 + a 3 + 4a 2 c-db 2 -bad + d 2 c 
+ ad 2 — bd 2 + ae 2 — dbe + 2ead + 6dm 2 a 2 ) 
( —a m — 2a — d + an + d) C m +2 (n) = 

with respect to m, with initial values C n (n) = l,C n+ \(n) = 0. Furthermore the recurrence 
equation 

(d + 2 an + 2 a) (n + 2)(-n 2 ab 2 - d 2 bn - dbe + 2a 2 n 2 e + 4n 2 a 2 c 

+ 2dn 3 a 2 + d 2 an 2 - db 2 n + a 3 n 4 + d 2 c + ae 2 - dan 2 b + 2dane + 4dcna) 
( n + 1) (an — a + c? + a m ) (—am — d + an + d) C m (n) — (d + 2an) 
(d — a + 2 an) (d + a + 2 an) (n + 2 )(— 2 ead — 2na 3 m 2 + edd 
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o cy — O O — Q O O 

+ 2 a nm + amd — am d +d nd — 2ea m + an dd — la m n 

— 2 a 2 mdn 2 — d 2 bn — a 2 m 2 d — an db + 2 a 2 n 2 e + 2 a 3 mn 2 + 3 a 2 n d 
+ 7a 2 n 2 d + andd + 2ead — dmd 2 + 4dn 3 a 2 — 4a 2 en 2 

-ed 2 + 2 a 2 ne + d 2 d + 2d 2 an 2 + 2dbd+2a 3 n 4 + ad 2 - bd 2 

— dan 2 b + a 2 md + 2dane + 4a 3 n 3 + 2a 3 n 2 — 2dmand — 2 a 2 m 2 n d 
+ 2 a 2 e m 2 + 2 a e m d — 4 ae n d — 2 ae d + 2 a 2 n dm 

— 2a 2 mdn — amdd — dbmd — a m 2 bd + ambd — 4ea 2 n 
+ 2dbnd + 2an 2 bd + 2anbd + 3and 2 )C m (n + 1) + 

an + 2 a) (d — a + 2 a n) (d + a + 2 an) ( —m + n + 2 ) 
(d + am + a + an) C m (n + 2) = 
with respect to n is valid. 

Next let cf(x) + f(x) = <j{x) + t(x), hence a = a, b = b + f,c = c + g, d = d — f,e = e — g for 
some constants /, g. Then the relation (22) is valid, C m (n) satisfying the recurrence equation 

( — d + f — 2 a m ) ( — d + f — 2 am — 2 a) 2 ( — d + f — a — 2 am) ( —d + f — 3a — 2 am) 
( — m + n) (an — a + d + am) C m (n) — ( — d + f — 2am — 2a) ( —d + f — a — 2am) 
(-d + f - 3a - 2am) (m + l)(2ea 2 m - 2a 3 m 2 n - d 3 + 2a 2 gm + 2ead - ad 2 

— bd 2 + d 2 bn + d 2 an 2 + a 2 n 2 d - 2a 2 n 2 e + 2a 2 ne - a 2 nd + dan 2 b - andb 

— 2dane — 2aef — dan f + 2am n +2anmd —2a m —am bd 

— ambd — la m d — 3a md — 3am d — 4amd +/ m — 4a m —2a m 

0-00 'i 

— amfb — 2a mfn +2a n dm + 2a fmn — 2a mn + 2a n m — md 

'i o 

+ / ban + 2dg an + 2a m nd+d n + 2a em +a nf — dbnf — d nf 

— a 2 n 2 f — an 2 b f — 2 g a 2 n + 2 g a 2 n 2 — 2 m f an d + m f d 2 + 3 / a d — 2 f 2 a 
-2f 2 d + 2fd 2 + 4am 2 fd + 2deam + d 2 g + 5a 2 mf + 9a 2 m 2 f-md 2 b 

— f e d — f g d — m f 2 d + 8 am f d + f 2 e — 2 a em f + 2 a d g m — 2 a f g m 

— 4a 2 dm 3 + 4 a 2 f m 3 + dfb + f 2 bm-3am 2 f 2 + 2a 2 m 2 g - am 2 f b - 6 f 2 am 
+ f 3 ) C m+ i(n) - ( -d + f - 2am) (m + 1 )(4ea 2 m + 8a 2 cm - 2b 2 am 

+ 4a 2 gm + 2ead + 4dca-db 2 + d 2 c + ad 2 -bd 2 -b 2 a + ae 2 + 2a 2 e 
+ 2 a 2 d + 4 a 2 c + a 3 - dbe - bad - 2a ef + 6 a 3 m 2 + 4a 3 m- am 2 bd 

— 2 a mb d + 6 a 2 m 2 d + 6 a 2 m d + a m 2 d 2 + 2 a m d 2 + 4 a 3 m 3 + a 3 m 4 
-b 2 am 2 +4a 2 cm 2 -2amfb + 2a 2 em 2 -mfd 2 -3fad + f 2 a + f 2 d- fd 2 

— 3a m 2 f d + 2 d e am + d 2 g — 6 a 2 m f — 6 a 2 m 2 f — md 2 b + 2dg a — fed 

— 2 f g a — fgd + mf 2 d — 6 am f d + f 2 e — 2aemf + ag 2 +4adcm + 2adgm 

— 4 a f cm — 2 a f g m + f 2 c — 2 d f c — 2 a e g + db g + f b e — f b g + 2 a 2 d m 3 

— 2 a 2 / m 3 - db 2 m + / b 2 m + f 2 bm + am 2 f 2 + 2 a 2 m 2 g + f 2 b - am 2 f b - a f b 
+ 2f 2 am-2a 2 f + fb 2 + 2a 2 g-4afc)(m + 2)(-d + f-am-an-a) 
(—am — 2a + f + an) C m+ 2 (n) = 

with respect to m, with initial values C n (n) = 1, C n+ \(n) = 0. Furthermore the recurrence 
equation 

(n + 2) (d + 2 an + 2a)(2dn 3 a 2 + a 3 n 4 + d 2 c + a e 2 - db e - db 2 n - d 2 bn 
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— n 2 ab 2 + d 2 an 2 + 2 a 2 n 2 e + 4n 2 a 2 c — dan 2 b + 4dcna + 2dane) 

(n + 1) (an — a + d + am) ( —a m + f + an) C m (n) — (d — a + 2an)(d + a + 2an) 
(d + 2an) (n + 2)(-2ea 2 m - 2 a 3 m 2 n + d 3 + 4dn 3 a 2 - 2 e a d + 2 a 3 n 4 
+ ad 2 + 6d 2 + 4a 3 n 3 + 2a 3 n 2 + d 2 6n + 3d 2 an 2 + 4and 2 + la 2 n 2 d 

— 2 a 2 n 2 e — 2a 2 ne + 3a 2 nd+dan 2 b + andb — 2dane + 2ae f — dan 2 f 

Q O O O 00 

— 2a m n —2anmd —am bd + ambd — a m d + a md — am d 

+ 2 a 2 m f n 2 - 2 a 2 n 2 d m + 2 a 2 / m n + 2 a 3 m n + 2 a 3 n 2 m - m d 3 - 2 / ban 

00 Q 00 n On On 

— d f an + 4dg an — 2 a m nd + d n + 2a em — 2dbn f — d n f — 2 an b f 

o o o 

+ 4ga n + 4ga n + 2m / and + m / d + m f bd — f d + 2deam + d g 

— m d 2 b + 2 d g a — f e d + am f d — 2 a em f — 2 d f b) C m (n + \) — (d — a + 2 an) 
(d + a + 2 an) (d + 2 an) 2 (d + 2 an + 2 a) ( —m + n + 2 ) ( — d + f — am — an — a) 
C m (n + 2)=0 

with respect to n is valid. □ 

Note that the recurrence equation for (7(2;) = cr(x) given in Theorem 3 reduces to two terms, 
and hence can be represented by hypergeometric terms, for the connection between Charlier 
polynomials (P n (x) = Cn\x),Q m (x) = Cm\x)), between Meixner polynomials (P n (x) = 
mn ,lx \x),Q m (x) = mm ^(x)), and between Krawchouk polynomials (P n (x) = k^\x 7 N), 
Qm(x) = km\x, M)), We will consider these and more cases by another method in § 6. 

3 Hypergeometric Representations: Continuous Case 

Note that by P^\x), C%{x), L { n a) {x), H n {x), B { n a \x) we denote the monic Jacobi, Gegen- 
hauer, Laguerre, Hermite and Bessel polynomials. Their non-monic counterparts have the 
standardizations (see [1], (22.3), and [2]; Al-Salam denotes the Bessel polynomials by Yn a \x j) 



system 


P^\x) 


CZ(x) L ( n a \x) H n (x) Bi a) (x) 




l / 2n + a + $ \ 
2" I n ) 


(a) n 2 n (-1)" on (n+a+l) n 



We get 

Theorem 4 Let P n (x) be a monic polynomial system given by the differential equation (2) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e. Then the power series coefficients C m (n) given 
by 

n 

P n (x) = C m (n) x m (28) 

satisfy the recurrence equation 

(m — n)(an + d—a + am)C m (n) + (m+l)(bm + e)C m+ i(n)+c(m + l)(m + 2)C m+ 2{n) = . (29) 
In particular, if c = 0, then the recurrence equation 

(m — n)(an + d — a + am)C m (n) + (m + l)(bm + e)C m +i(n) = 
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is valid, and we have the hypergeometric representation 



P n {x) 



(t)» OH 



valid for a / 0, or 

Pn(^) = 

valid for a = 0, b ^ 0, or finally 

PM = {I 



— 1 -2^11 



"I I 1*1 



-n 



-n 



d+(n— 1)« 



b\ 



a 



1*0 



-n 



(30) 



(31) 



(32) 



valid for a = 0, b = 0. 

Therefore, the classical continuous orthogonal polynomials and their monic counterparts have 
the following hypergeometric power series representations: 

-n, n + a + (3 + 1 
a + 1 

^ 2n + a + /3 j ^ ~ T 



p(^) ( ,) = ( n ; a ) 2 Fi( 



1 — a; 



-n, — n — a 



-2n-a-/3 



1-x 



/ 2n + a + /? ^ /'z + 1 
V n 



-n, n + a + /3 + 1 
/3+1 



1 + a; 



n / -n,-n-/3 
2*1 ( 



-2n -a- (3 



1+x 



C%(x) = x n 2 F 1 \ 



-n/2,-n/2 + 1/2 
— n — a + 1 



1 



9 I ' 



(«)»2^ / -n/2, -n/2 + 1/2 
C rc W = 2^1 



L^(x) = (l + a) n (-l) n 1 F 1 
LP(x) = ( n l a ) 1 F 1 ( i 
H n (x) = x n 2 F 



—n 
1 + a 



-n — a + 1 



x \ = x n 2 F 



9 I ' 



-n, —n — a 



-n 
+ a 



x 



-x) n / — n, — n — a 

— j 2*0 

n! \ — 



a; 



2-^0 



4 Q) (*) 



n/2,-n/2 + 1/2 
-n/2, -n/2 + 1/2 
n, n + a + 1 



2-^0 



4^) = 2*0 



(n + a + l) n y 
n, n + a + 1 



z"l*l 



a; 
2 



(n + a + 1) 
2" 



— n 


2 


—2n — a 


X 


—n 


2 


—2n — a 


X 



These results are all particular cases of the recurrence equation (29). 
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Proof: Substituting the power series (28) into the differential equation, and equating the 
coefficients yields the recurrence equation (29). 

For c = this recurrence equation degenerates to a two-term recurrence equation, and hence 
establishes the hypergeometric representations (30)-(32), using the initial value C n (n) = 1. 
A shift in the rr-variable then generates the representations for the Jacobi polynomials. The 
two points of development x\ = 1 and X2 = —1 correspond to the zeros of a(x). 
Note that some of the hypergeometric representations correspond to each other by changing 
the direction of summation. 

The other representations follow by substituting the particular parameters a, 6, c, d, and e into 
the recurrence equation (29), and using the initial value C n (n) = k n (or C n (n) = 1 in the 
monic case). □ 

We would like to mention that the recurrence equation (29) carries the complete information 
about the hypergeometric representations given in the theorem. 

The method described results in four different hypergeometric representations for the Jacobi 
polynomials. Many more hypergeometric representations exist, but the algorithmic procedure 
presented finds power series representations only. For example, the representation (see e.g. [1] 
(22.5.45)) 



^ (x) = (-^)(£zl)' 



/ —n, —n — a 

2Fl[ 



p + 1 



X + l 



cannot be discovered by this method. 

The method was able to find hypergeometric series representations with point of development 
xq = for the Gegenbauer polynomials which are specific Jacobi polynomials, but failed in 
the Jacobi case, though. One might ask whether such a representation exists. This question 
can be completely answered by an algorithm of Petkovsek [18]. Petkovsek's algorithm finds 
all hypergeometric term solutions of holonomic recurrence equations, i.e., homogeneous linear 
recurrence equations with polynomial coefficients. Using the recurrence equation (29), an 
application of Petkovsek's algorithm proves that the Jacobi polynomials do not generally have 
a hypergeometric series representation at the origin. 

Note that the method of the last section, although more complicated, does also give the 
recurrence equation (29), and hence the above results. 



4 Power Representations 

Whereas in the last section we considered the specific connection coefficient problem for 
Qm(x) = x m , in this section the opposite problem, having P n (x) = x n , is studied. 
In many applications, one wants to develop a given polynomial in terms of a given orthogonal 
polynomial system. In this case handy formulas for the powers x n are very welcome. 

Theorem 5 Let Q m (x) be a monic polynomial system given by the differential equation (2) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e. Then the coefficients C m (n) of the power 
representations 

n 

x n = C m (n)Q m (x) 

satisfy the recurrence equation 

(n — m) (d + 2am) (d + 3a + 2am) (d + a + 2am) (d + 2am + 2a) 2 C m (n) 
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+ (de + bd + 2dbm + 2a m 2 b + 2amb + 2ean — dbn) (d + 2 am + 2 a) 

(m + 1) (d + 3a + 2am) (d + a + 2am) C m+1 (n) - (rn + 2)(-4a 2 cm 2 (33) 

+ a b 2 m 2 + 2ab 2 m — 4acmd — 8a 2 cm + mb 2 d — ae 2 —d 2 c + bed — 4a 2 c 

— Aacd + at> 2 + t> 2 d)(am + an + a + d)(m + l)(d + 2am) C m+ 2(n) = . 

If c = 0, then the recurrence equation 

(n—m) (d+2am) (d+a+2am) C m (n) + (m+1 ) (bm+e) (am+na+d) C m+ i(n) = (34) 

is valid, and we get (a b / 0) 

Therefore, the following representations for the powers in terms of the classical continuous 
orthogonal polynomials are valid: 

(1 - xY = 2" T(a + n + l)± + ^ + ^ + l)T(a + + m + 1) p{ p) 

V ' y ' ^ T(a + m+l)T(a + f3 + n + m + 2) V ' rn y ' 

(see e.g. [20], 136, Eq. (2), or [17], § 5.2.4; note the essential misprint in this formula!), 

(i + ,)» = 2» ros + n + 1) t (-1)- '^^n'^tft"!!! <->- • 

1 (p + m + 1) 1 (a + p + n + m + 2) 

T n = L V (-»/2) fc (-n/2 + l/2) fc (-n-«) fc ~ Q 

^ (-n/2 - a/2) fc (-n/2 - a/2 + l/2) fc k\ \ a) ^-2kK ) , 



x 



nl L v^ J {-n/2 - a/2 + l) k (-n - a) k . , 
" («)»2» fc tj (-n/2 - a/2) fc *! 1 ij G "" 2 ^ j 



_ n! L ^ J n + a-2fc 
(see e.g. [20], 144, Eq. (36), or [17], § 5.3.4), 

s B = (l + a) B £J=^ ^ (; + «)(-!)•» LW (a; ) 

m=0 v ~r M to= q 
(see e.g. [20], 118, Eq. (2), or [17], § 5.5.4), 

_ L ^ J (-n/2), (~n/2 + 1/2), ~ 
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r n _ ^ (-n/2) fc (-n/2 + 1/2), _ n! L ^ J 1 



A:=0 ' fc=0 

(see e.g. [20], 110, Eq. (4), or [17], § 5.6.4) 

x n (-2)" " 

(a + 2)„ ^ (n + 2 + «) m m! 



(~")m (a/2 + l)m (a/2 + 3/2) m ~ (a) 



("2) n A (-n) m ( a +l) m (a/2 + 3/2) TO (a) 



(« + 2) n ^ (n + 2 + a) m (a/2 + l/2) m m! 
- f-2) n V f 2m + a + 1) T(ft + m + 1) (a) 

m=U v ' 

(see [2], (7.5); note the essential misprint in this formula!; compare [20], 150, Eq. (7)). 

Proof: In § 2 it was shown how one obtains three essentially different cross rules for the 
connection coefficients between P n (x) and Q m (x). We modify this method here. For Q m (x), 
we have the differential equation 

a{x) Q" m (x) + f{x) Q' m (x) + A m Q m (x) = 

with a(x) = ax 2 + bx + c, and the derivative rule 

a{x) Q' m (x) = a m Q m +\ (x) + (3 m Q m {x) + 7m Qm-i {x) , 

and it is easily seen that our current P n (x) = x n satisfies any of the derivative rules 

a(x) P' n (x) = anP n+ i(x) + bnP n (x) + cnP n -i(x) . (36) 

Hence in our situation, we get the two cross rules (23) with a n = 1, b n = c n = 

C m {n + 1) = a m -i C m _i(n) + b m C m {n) +c m+ i C m+1 {n) (37) 

and (25) with a n = l/(n + 1), b n = c n = 

C m (n + 1) = 2 m _i C m _i(n) + 6 m C m (n) + c m+ i C m+ i(n) (38) 



n + l 

which we had deduced in § 2. Using the derivative rule (36), we obtain the third cross rule 

a n C m (n + 1 ) + b n C m (n) + c n C m (n - 1 ) = a TO _ i C m _i (n) + /3 TO C m (n) +7 m +i Cm+i («) • (39) 

To receive the recurrence equation (33), we use Theorem 1 writing the cross rules in terms 
of a, b, c, d, and e, only. Then by linear algebra we eliminate the variables C m (n + 1) and 
C m (n— 1) to obtain a pure recurrence equation w.r.t. m. (Similarly by elimination of the 
variables C m _i(n) and C m +i(n) a pure recurrence equation w.r.t. n is obtained.) A shift by 
one gives (33). 

If c = 0, then the recurrence equation has still three terms, unfortunately. But since for c = 
in neither of the three cross rules (37)-(39) the variable C m (n — 1) does occur, we can do a 
similar elimination, this time eliminating the variables C m (n+1) and C m _i(n), leading to the 
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first order recurrence equation (34). Hence the hypergeometric representation (35) follows. 
The power representations for the Jacobi, Laguerre and Bessel polynomials are special cases 
thereof. 

In the case of Hermite and Gegenbauer polynomials, (33) contains only the two terms C m (n) 
and C m+ 2(n), which leads to the desired representations. □ 

Note that, again, the recurrence equation (34) carries the complete information about the 
hypergeometric representations given in the theorem. 

As an immediate consequence of the above theorem, we get the following connection coefficient 
results. 

Corollary 2 The following connection relations between the classical orthogonal polynomials 
are valid: 

P^(x) - V(2m + 7 + /3 + l) r(n + /3 + 1) nn + m + a + (5 + l) 



m=0 



r(m + 7 + 8 + 1) (a- j) n - m (7:/3) 

' - r m \ x ) 



T(n + m + j + 8 + 2) (n - m)\ 

(see e.g. [4], (13)), 

P^W = ±{-lY- m {2rn + a + 5+l) ^ n + a + l \ nn + m + a+ /\ 1) ■ 
n ' ' V 'T(m + a+l) T(n + a + 8+1) 



l ,,,„ , ... , 1) Tin + a + 8 + 1) 

T(m + a + 5 + 1) {8- <J)„_ r 



T(n + m + a + S + 2) (n — rn)\ 
(see e.g. [6], (2.8)), 

ra , , _ T(/3) L ^ J ( n -2k + 8)T(k + a-8)T(n-k + a) p 
n[) ~T(a)T(a-8)^ k\ T(n - k + 8 + 1) Vm1 J 

(see e.g. [5], (3.42)), 

4" ) w=E ( r /3) v L -^ 

(n — ml! 

(see e.g. [20], 119, Eq. (2)), 

„(«), , (-l) n (a-B) n " (-n) m Q9+l) m Q3/2+3/2) m (n+a+l) m (/3) 

" M " W+2) n ^ (n+2+8) m (8/2 + l/2) m (8-a+l-n) m m\ { ' m{> 

- f f l) m (2m+8+l) (- n )mnP+m+l)(n+a+l) m T(P-a+l) (0) 

~ io ^ m+ ^ +1 ^ m! r(n + m + /3 + 2)r(m-n + /3-a + l) m 1 j ' 

(see [2], (8.2); note the essential misprint in this formula!). 

Proof: We want to find the coefficients C m (n) in the relation (22) 

n 

P n(x) = ^2 C m{n)Q m {x) . 
m=0 
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Combining 

P n {x) = 'Y^A j (n) x 3 and x 3 = E B m (j) Q m (x) 

iez mez 



yields the representation 



P n(x) = E E A i(") B ™U) Qrn(x) 



and interchanging the order of summation gives 



C m (n) = E^jW^mO') • 

Similarly, if (as in the Gegenbauer case) 

P n {x) = E Aj(n) a;"" 2i and ^ = E B m(j) Qj-2m(x) 

jez mGZ 

one gets 

= E 4/( n ) ^™-i( n - 2 i) 

iez 



with 



P„(a;) = E D m{n)Q n -2m{x) 

m=0 



Since the summand F(j, m, n) := j4j(n)S m (j) turns out to be a hypergeometric term with 
respect to (j, m,n), i.e., the term ratios F(j + 1, m, n)/F(j, m, n), F(j,m + l,n)/F(j,m,n), 
and F(j,m,n + m, n) are rational functions, Zeilberger's algorithm ([26], [11], see e.g. 

[9]) applies and finds recurrence equations for C m (n) with respect to m and n. 
In all cases considered, Zeilberger's algorithm finds recurrence equations of first order with 
respect to m (as well as for n). The given representations follow then from the initial value 
C n (n) = k n /k n . □ 

For some applications, it is important to know the rate of change in the direction of the 
parameters of the orthogonal systems, given in terms of the system itself. By a limiting process, 
these parameter derivative representations can be obtained from the results of Corollary 2. 

Corollary 3 The following representations for the parameter derivatives of the classical or- 
thogonal polynomials are valid: 

an— l 

da ±n ~ ^a + 6+l+m + n 

a + (3 + 1 + 2m ((3 + m + !)„_ 
n — m (a + f3 + m + 1) 

(see [7], Theorem 3), 

/ 2m + a + (3 \ 

1_ p(a,p )(x) = y [ m I « + /3 + l + 2m (/3 + m + l) n - m p(a , p) 

da n £^ Q n-m I 2n + a + f3 \ a + /3 + 1 + m + n {a + /3 + m + l)„_ m m 

\ n ) 



pW\x) = E 5— r (Pt p) (x) + 

-P^(z)) 
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n-1 



= E „ +l , + ; +m+ „ -(^M + 

3 + 1 + 2m (a + m + 
n — m (a + f3 + m + l) n _ 



(see [7], Theorem 3), 



i^n-m « + 1 + 2m (« + m + l)»_ m p ( Q ,/3)/ A 



/ 2m + a + /? \ 

lpW) w = (~2)"- m I m J a + /3+l + 2m (a + m + l) n - m p {a ,p )(x) 
9(3 n ^ n - m (2n + a + {3\a + (3 + l + m + n{a + (3 + m + l) n _ m m 1 J 



±c°ix) = r( 2(1+m) + - \c°ix) 

da n{ ] ^ \{2a + m){2a + l + 2m) 2a + m + n) n{ 1 



n ^ 2(1 + (-!)"-"» )(a + m) 
+ (2a + m + n)(n-m) 



(see [12], Theorem 10), 



JLr a (r\ - ^™ n! (i + (-i)"- m )(ft + m) ~ 

da " ^ (a) n m! (2a + m + n) (n - m) m{ ] 

= V " ! »-2fc + « ~ 

jfo (a + n- 2k) 2k 4 k (n - 2k)\ k (n - k + a) n ~ 2M ' ' 

o ra— 1 ^ 

oa '„ n — m 



m=0 



(see [12], Theorem 10), 



da " 1 j ^ n-m m! m 1 j 

m=0 



o n— 1 i 

^-K(^) = E 7 • fe*) + 

da n ' ^a + n + m + 1 \ ny ' 

m=U 



( _ ir . jS±£+l " ! asw) . 

(n — m) (a + m + l) n -m m - 1 

_£_ B a ( x ) = V — - - £ a (a;) 

da " ^ (n - m) (a + n + m + 1) (a + 2m + 2) 2 „_ 2m _i m! mV 

Proof: Given the connection relation 

n 
m=0 
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we build the difference quotient 



P%{x)-Pj{x) " C m (n;a,0) Pj(x) 

a-0 ^ a-0 m[X) a-0 

n-1 



C n (n;a,P) - 1 p ^ C m (n;a,P) p 



so that with (3 — > a 



W = Jim - 1 CW + £ J im K(x) ,40) 

oa p-ya a — p ^-r, p-ya a — p 

1 m=0 ' 

since the systems P®{x) are continuous w.r.t. a. This gives the results. □ 

Note that for monic polynomials (and moreover if k n does not depend on a as in the Laguerre 
case) the first limit in (40) equals zero. Hence the parameter derivative representations are 
simplest in such a case. 

5 Hypergeometric Representations: Discrete Case 

By h ( n'^(x,N) and Q n (x; a, (3, N) we denote two commonly used standardizations of the 
Hahn polynomials (see [19], and [23]), and by mi 1 '' 1 '(i), kn\x,N) and c^\x) the Meixner, 
Krawchouk and Charlier polynomials are denoted, respectively. They have the standardiza- 
tions 



system 



h { n a ' P) (x,N) Q n (x;a,(3,N) m^\x) h£\x,N) <&\x) K^\x) 



k n 



(a + (3 + 2n\ {a+ p +n+ i )n I _ A" ± /A" 

V n ) {-N) n { a +i) n \ L n) n\ \ n) a 



The polynomials t n (x,N) := hn'°\x, N) are the discrete Chehyshev polynomials. The poly- 
nomials K^C ,/3 \x) given in (21), are not orthogonal, but satisfy the difference equation 

AVy(x) + {ox + (3) Ay(x) + X n y(x) = . 

The monic counterparts of the discrete systems will be denoted by hi?'^ (x, N), Q n (x; a, (3, N), 
t n {x,N), fhti'^ix), kli \x,N) and c^{x), respectively. Observe that therefore by h n we do 
not denote the Hahn-Eberlein polynomials h^ ,v \x,N) as in [19]. 

In the continuous case, we looked for power series representations, i.e., we set Q m (x) = x m . 
The corresponding choice in the discrete case is a representation in terms of the falling factorial 

Q m {x) =x^:=x(x-l)---(x-m + l) = (x-m + l) m = (-l) m (-x) m . 

We get 

Theorem 6 Let P n (x) be a monic polynomial system given by the difference equation (4) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e. Then the series coefficients C m (n) given by 

n 

P n (x) = C m (n)x& (41) 

m=0 
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satisfy the recurrence equation 



(an + am — a + d) (n — m) C m (n) 

+ (m + 1) (an 2 — 2a m 2 — an — am + nd — 2 dm — bm — d — e) C m+ i (n) 

— ( rn + 1 ) ( rn + 2 ) ( a m 2 + 2am + dm + bm + a + d + b + c + e) C m+2 (n) = 

If c = 0, then the recurrence equation 

(n — m) (am + d + an — a) C m (n) — (m + 1) (a m 2 + mb + md + e) C m+ \ (n) = 

is valid, and we have the hypergeometric representation 

(d_A ( b+d+^{b+df-Aae \ ( b+d-^ \b+d) 2 -4ae " 



Pn(x) 



( A\ ( d - a \ 
\2a) n \ 2a ), 



4, 



— n, —x, n — 1 + - 

'3^2 I b+d+^J(b-d) 2 -4:ae b+d-^ \b-d) 2 -4ae 



2a 



2a 



valid for a / 0, or 



Pn(x) 



,b + dj n 

valid for a = 0, b + d ^ 0, or finally 



1 + 



-n, —a; 

e 

6+d 



6 + d 



2-fb 



-n, —a; 



valid for a = 0, d = —b. 

Therefore, the classical discrete orthogonal polynomials and their monic counterparts 
the following hypergeometric series representations: 



hW\x,N) 



l) n , / -n, -x,n + 1 + a + 8 

J-(8 + l) n (N-n) n3 F 2 [ 



(3 + 1,1 — N 



1 



(see e.g. [19], p. 54, Table 2.4), 

*£■*><,,»)- < 1 + «-< 1 - JV > 



(1 + n + a + 3) n 



( — n, — x, n + 1 + a + 3 

3^2 

1 P+l,l-N 



t n (x,N) = {-l) n {N - n) n ■ 3 F 2 \ 



(l/2)„ 4« 



Q n {x;a,P,N) = 3 F 2 



—n, —x, n + 1 
1,1 -iV 

■n, —a;, n + 1 
1,1 -JV 



-n, — rr, n + 1 + a + /? 
a + 1, —N 
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(see e.g. [23], 1.5) 



Q n (x;a,(3,N)=h^ a \x,N + l) , 



2^11 



m^\x) = { 1 ) n2 F l [^ ^ X 
(see e.g. [19], p. 54, Table 2.4), 

m^(x) = ( 1 ) n 

V A* 

(see e.g. [19], p. 54, Table 2.4), 

(x,N) = (-N) n p n 2 F l (~ 



-n, — x 
7 



i-I 



-n, —a; 
-N 



c£>{x) 



-n, —x 



2-f o 



n, —a; 
—N 

-a 



(see e.g. [19], p. 54, Table 2.4), 

4")(x) 



-n, —a; 



These results are all particular cases of the recurrence equation (43). 

Proof: Substituting the series (41) into the difference equation, and equating the coefficients 
of the falling factorials yields the recurrence equation (42) which had been obtained by Lesky 
[16]- 

This conversion can be easily done using a computer algebra system by bringing the given 
difference equation into the form (5), expanding it, and replacing any occurrence of Ay(x) by 
(m + 1) C m +i, any occurrence of a product xy(x) by C m _i + mC m and any occurrence of a 
shift y(x + 1) by C m + (m + 1) C TO+ i since 



xx m = x m±l + mx m 



and {x + \) m = x m + mx^- 



Iteratively for all nonnegative integers j, k any of the terms x^ A k y(x) and x^ y(x + k) can be 
replaced by these rules. Note that this method can also be applied for higher order difference 
equations with polynomials coefficients. 

Different from the continuous case, the recurrence equation (42) does not degenerate to a 
two-term recurrence equation for c = 0. To get (43), nevertheless, we must use a different 
approach. One possibility is to apply Petkovsek's algorithm to the recurrence equation (42), 
leading to (43). 

Another possibility is to modify the method which will be used in the next section to deduce 
representations of the falling factorials in terms of discrete orthogonal systems. This method 
yields (43) directly. 
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As soon as (43) is deduced, the initial value C n (n) = 1 gives the hypergeometric representations 
(44)-(46) which include all other representations by substituting the particular parameters 
a, b, c, d, and e. □ 

We would like to mention that, again, a single recurrence equation, (42), carries the complete 
information about the hypergeometric representations given in the theorem. 
Note furthermore, that the radicals in (44) do only occur by the representation used: the 
radical factors come in pairs whose product is radical-free. Note that the computation which 
gives (44), answers a question raised by Koornwinder [15]. For more examples of this type see 

[13]- ' . , , 

Our method was able to find hypergeometric series representations for the particular case 
c = 0. This is the most important situation since all the classical discrete orthogonal families 
are of this type, corresponding to the fact that their discrete support has zero as left boundary 
point (see e.g. [19], Tables 2.1-2.3). 

By construction, all the series representations determined have an upper parameter —x. The 
question remains, however, whether or not such a hypergeometric series representation might 
be valid for c / 0, too. In general, the answer is no. Petkovsek's algorithm shows that the 
recurrence equation (42) does not generally have a hypergeometric term solution. 
Note that the hypergeometric representation (21) for K^ ,/3 \x) is not of this type, and cannot 
be obtained by the given method. By Petkovsek's algorithm there is no representation (41) 
with a hypergeometric term C m (n) for these polynomials. 

6 Falling Factorial Representations 

Whereas in the last section we considered the specific connection coefficient problem for 
Qm(x) = x—, in this section the opposite problem, having P n (x) = x—, is studied. 

Theorem 7 Let Q m (x) be a monic polynomial system given by the difference equation (4) 
with a(x) = ax 2 + bx + c, and t(x) = dx + e. Then the coefficients C m (n) of the falling 
factorial representations 

n 

E C m (n)Q m (x) (47) 

m=0 

satisfy the recurrence equation 

(2ma + a + d) (2ma + 3a + d) (2ma + 2a + d) 2 (2ma + d) (n — m) C m (n) + 

( 2 m a + a + d ) (2ma + 3 a + d) (2ma + 2a + d) ( m + 1 ) (2 m 2 n a 2 — 2 m 2 a 2 

+ m 2 ad + 2 m 2 ab + 2mna 2 + 2mnad — 2m a 2 — mad + 2m ab + md 2 

+ 2 m db + na d + 2nae — ndb — ad + db + de) C m+ \(n) 

+(m + 1 ) (2ma + d)(m 4 a 3 +4m 3 o 3 + 2rn 3 a 2 d + 6rn 2 a 3 + 6m 2 a 2 d 

o o o o o — 2 o o — 2 Q o — 

+ 4m a c + 2m a e + m ad — m adb — m ab + 4 m a + 6 m a d 

_2 _ —2— — — 2 — _ — _ _ — _ 

+ 8ma c + 4ma e + 2mad — 2madb + 4madc + 2made 

— 2m at> 2 — md 2 ~b — mdt> 2 +a 3 + 2 a 2 d + 4 a 2 c + 2 a 2 e + a d 2 

— a db + 4a dc + 2 a de — at> 2 + a e 2 — d 2 b + df 'c 

-db 2 -dbe)(m + 2) (ma + na + a + d) C m+2 (n) = . (48) 
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If c = 0, then 

= (d + a + 2am) (d + 2am) (-n + m)C m {n) (49) 
— (an + d +am) (m + 1) (am 2 + md + mb + e) C m +i(n) . 

Therefore, the following representations for the falling factorials in terms of the classical dis- 
crete orthogonal polynomials are valid: 

n= A ((3 + l) n (l-N) n (-ir(l + a + (3 + 2m)(-n) m (l + a + (3) m 

^ (a + (3 + 2) n (l + a + (3)(n + 2 + a + (3) m ((3 + l) m (l-N) m m { ' ; ' 

n = " Q3 + 1)„ (1 - N) n (-1)" (-n) m (a/2 + g/2 + l) m (q/2 + (3/2 + 3/2) m 4™ ~ (a ^ 

^ (« + /3 + 2)„ (n + 2 + a + / 3) m (/3 + l) m (l-7V) m m! m ^' ; 



m=0 



„ _ A (l + a)n(-N) n (-ir (a + g + 1 + 2m) (-n) m (1 + a + 
* "^0 (a + /3 + 2)„ (« + /3 + l) (n + 2 + a + /3) m m! ^'"'^ 

(compare [8], (4.2)-(4.3)), 

(l-N) n (-l) n " (_ n ) m (i + 2m) 



(1-AQ W (-1)" A (-n) m (3/2) m 4™ ~ 

n + 1 ^ (n + 2) m (l-iV) m m! tm ^ 7Vj 

" (-l) n (7)nfe) n (-«) m , , 
^0 (7) m "»! m {U 

" ("If (7)n (t^)""" 1 (-n) m . . 

(-N) m m\ 
m! 



*= E ( " 1)w( " /i) 7" i ( " n)m ^(x) . 

m=0 m ' 



Proof: In § 2 it was shown how one obtains three essentially different cross rules for the 
connection coefficients between P n (x) and Q m (x). We modify this method here. For Q m (x), 
we have the difference equation 

a(x) AVQ m (x) + t(x) AQ m (x) + A m Q m (x) = 
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with <j[x) 



ax 2 + bx + c, and the difference rule (11) 



(a(x)+T(x))AQ Qm+l{x) + (fi 



TO 



•^m) Qm(^) "I - 7m Qra-l(^) j 



and it is easily seen that our current P n (x) = x— satisfies any of the difference rules 




(50) 



Hence in our situation, we get the two cross rules (23) with a n = 1, b n = n, c n = 



C m (n + 1) = a m -\ C m _i(n) + b m C m (n) + c m+ \ C m+ i(n) 



(51) 



and (25) with a n = l/(n + 1), b, 



'n 



= 



1 



a n C m (n + 1) 



a m -i C m _i(n) + b m C m (n) + c m+1 C m+ i(n) 



(52) 



n + 1 



which we had deduced in § 2. Using the difference rule (50), we obtain the third cross rule 



To receive the recurrence equation (48), we use Theorem 1 writing the cross rules in terms 
of a, 6, c, d, and e, only. Then by linear algebra we eliminate the variables C m (n + 1) and 
C m (n— 1) to obtain a pure recurrence equation w.r.t. m. (Similarly by elimination of the 
variables C m _i(n) and C TO+ i(n) a pure recurrence equation w.r.t. n is obtained.) A shift by 
one gives (48). 

If c = 0, then the recurrence equation has still three terms, unfortunately. For c = 0, we find 
a fourth cross rule to eliminate one more variable in the following way. Since the (second) 
difference rule (10) 

o{x) VQ m (x) = a m Q m+1 (x) + (3 m Q m {x) + 7 y m Q m -i(x) 

is valid, we can use the fact that for c = any of the difference rules 



and therefore the fourth cross rule 

anC m (n+ 1) + n (an + b) C m (n) = a m _ 1 C m _i(n) + (3 m C m (n) + 7 m+1 C m+ i(n) (54) 



Eliminating the variables C m (n+1), C m _i(n), and C m -\(n) from the four cross rules (51)-(54) 
gives the first order recurrence equation (49). This leads to the desired representations. □ 

Whereas we admit that all the shifted factorial representations of the theorem essentially 
were known ([8], [22]), our presentation unifies this development. In particular, the recurrence 
equation (49) carries the complete information about the falling factorial representations given 
in the theorem. 

Petkovsek's algorithm proves, again, that for the family K^C'^ (x) there is no representation 
(47) with a hypergeometric term C m (n). 

As an immediate consequence of Theorem 7, we get the following connection coefficient results. 



a n 




a(x) VP n (x) = anP n+ i(x) + n(an + b) P n (x) , 



is valid. 
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-h^Hx,N) 



Corollary 4 The following connection relations between the classical discrete orthogonal 
polynomials are valid: 

(2 + a + 5) n n\ 

(a + 5+l + 2m) (~n) m (1 + a + S) m (n + 1 + a + (3) m ,( a s) ( N) 

{a + 8 + 1) (l- N ) m (a+i) m (a+2+n+5) m (-n-(3+5+l) m ™ (X ' > ' [ ' 

r W)fT m = (« + l)n (1 ~ jVQn (1 + a + ff)n (I ~ S) n 

n [ ' j ^ (2 + « + 8) n (a/2 + (5/2 + l/2) n (a/2 + /3/2 + l) n 4" 

(-n) m (n + 1 + a + (3) m (a/2 + 5/2 + l) m (a/2 + J/2 + 3/2) m 4 TO ~ (a g) 

(l-7V) m (a + l) m (a + 2 + n + 5) m (-n-/3 + 5 + l) m m! m ^ j ' 

K ( ' ] ~h (2 + /3 + 7),n! 

(/3 + 7 + 1 + 2m) (-n) m (1 + g + 7 ) m (n + 1 + a + /3) m (-l) m 
(0 + 7 + 1) (l-JV) m (/3 + l) m (/3 + 7 + n + 2) m ( 7 -a-n + l) 

~ h (*,P)< N) = V ^ + 1)n (1 ~ (1 + " + g)n (« - 7)n (~1)" 

n 1 ' J i (2 + (3 + 7 )n (a/2 + /3/2 + l/2)„ (a/2 + (3/2 + 1)„ 4" 

(-n) m (n + 1 + a + /3) m (/3/2 + 7 /2 + l) m (/3/2 + 7 /2 + 3/2) m (-4) m 

(/3 + 7 + n + 2) m (l-7V) m (/3 + l) m (7-a-n + l) m m! m ^ j ' 

/^< Q )(x,Af) = (56) 

L ^ J (-n/2), (-( n -l)/2), (q- 7 ), ((N-n)/2) k (N-n+l)/2) k (- ra - 7 -l/2) fc - (7|7) 

(1/4 - 7/2 - n/2), (-n + 1/2 - a), (-n/2 - 1/4 - 7 /2), fc! 4* - 2fcl ' J ' 

h M ( iV) = (a + l) B (a + l/2) n (2 7 + l) n 
" 1 ' (7 + l)„( 7 + l/2)„(2a + l)„ 1 ; 

*=o ("T " n /2)fc ("7/2 - n/2 - 1/4), (-„-« + 1/2), (-7 - n/2 + 1/2), A! rfl ' J ' 

n ^ « « AH («+V2)n(27+l)n 

Q„(a;; a, a, JV) = 7^~TTi ' ( 58 ) 

(7 + l/2)n (2a + l) n 

L V (-n/2), (-(n - l)/2), (a - 7), (3/4 - 7 /2 - n/2), (-7 - n - 1/2), 

jfo (-7-n/2),(-7/2-n/2-l/4),(-n-a+l/2),(- 7 -n/2 + l/2),A;! ^-2^,7,7, ), 

n U n, R N\ X- (/5-^)n(-l) n 



m=0 (2 + a + S) n 
(a + 5 + 1 + 2m) (-n) m (1 + a + S) m (n + 1 + a + /3) m (-1)' 



(a + 5 + 1) (a + 2 + n + S) m {1-/3 + 5- n) m ml 

(compare [8], (4.1), (4.5)), 

n ( a KT\ f (<* ~ 7)n (ff + l)n 
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Q m (rr;a, J, iV) 



(/3 + 7 + 1) (P + l) m (P + l + n + 2) m ( 1 -a-n + l) m m\ W ' U,J 



(compare [8], (4.1), (4.5)), 



(compare [8], (5.5)), 

n / ,, \ n—m 



m=0 

(compare [8], (5.4)), 



(7,,) (a:) - f ( ^ V"™ (7)n(-n) m (-ir -(,,„),, 



(,, iv) = £ (p - g)— » ( — (*> 

(compare [8], (5.11)), 

(x, jv) = e (? - g)- m ( " Ar) ;_ ( "" )m( ," ir (*> *0 > 

(compare [8], (5.12)), 

" ( ' j "^ (iV-M-n + l) m m! &m 

4^) = E ^> - ^r m 4^) 

(compare [8], (5.16)), 

n ( — \ 

= E (-!) m - »)" ' 

/fl , N n (- n ) ( «(l-n)-0+6 \ 
K ^) {x)= (P^l\ a n-m ^ y j — l^K^[x) 

/_ \ / a(l-n)-/3+d' \ 

(,) = (^i) E 1 r — - ^ ^ • 

V a n ' ml 
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Proof: Combining the representations 



Pn{x) = A i( n ) xL and xl = E B m{j) Qm{x) , 

jez mez 

and using Zeilberger's algorithm, the method of Corollary 2 yields the results. 

The connection relations for the polynomials K^^\x) cannot be obtained by this method. 

Here Theorem 2 leads straightforwardly to the result. □ 

Although besides (56)-(58) the connection results were essentially known ([8], [3]), our de- 
velopment gives a unified treatment of them and makes new results like (56)-(58) easily 
accessible. 

Note that some of the representations are rather complicated. We suggest the idea to use the 
notation p f q for the summand of p F q , i.e. 



p ( u PP er 
P \ lower 



( upper 



X'. k 



k = \ lower 

With this notation, (55) could be written in the standardized hypergeometric notation 

' (3-6, 1-N,a+1 



ht®(x,N)= 3 fA 



2 + a + 6 



1; n 



" / -n, l + a + 6,n + l + a + (3,a/2 + 6/2 + 3/2, 1 
' ^^\l-N,a+l,a/2+6/2 + l/2,a+2+n+6, -n-0+6+1 



m=0 



l;m ht d) (x,N) . 



Finally, we deduce the parameter derivatives for the classical discrete orthogonal polynomials. 

Corollary 5 The following representations for the parameter derivatives of the classical dis- 
crete orthogonal polynomials are valid: 



d 



n-l 



1 



- • (h,W\x,N) + 



m =o oc + (3 + m + n + 
(-l)"" m (a+/3+l+2m) (l-7V+m)„_ m (p+i+ m ) f 



ht P) (x,N)) 



d_ 

da 



n-l 



(n - m) (a + (3 + 1 + m)„_ m 
(a+/3+l+2m) (l-7V+m)„_ m {0+l+m) n - m n! ~ (ce /3) 
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—Q n {x;a,(3,N) 



(a+(3+m+n + l) (n — m) (a+/3+l + 2m) 2n _ 2m m! 
1 



n-l 
V , 

' Vo; + /3 + m + n + l a + m 

m=u ' 

(« + /3 + 1 + 2m) (/3 + 1 + Tn) n — m n\ 



(n - m) (a + 1 + m) n _ m (a + /3 + 1 + m) 



— — -) • (Q n {x;a,0,N) + 
m + 1 / V 

Q m {x; a, 



s/3 



^ a,/3 )(a:,JV) 



n-l 

£ 

m=0 



1 



a + /3 + m + n + 



-.(h^\x,N) + 



a + /3 + 1 + 2m (1 - TV + m)„_ m (a + 1 + m) f 



n — m 



{a + (3+ 1 + m) ? 



■^(^JV)) , 
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— h^)(r .\n=V «+/3+l + 2m (l-N+m) n _ m {a + l+m) n _ m n\ ~ {a0) 

d(3 n K ' ' ^ (a+(3+m+n+l)(n-m) («+/3+l+2m) 2 „_ 2m m! m { ' h 



d 



Q n (x;a,(3,N) = ^ — ■ (Q n (x;a,P,N) + 

ap _ a + p + m + n+ l V 



m=0 



(-l)"- m (a + /3 + l + 2m) n! 



n — m (a + /3 + 1 + m) n _ m m! 

(see e.g. [12], Theorem 9), 

' m=0 v ' 

(see [12], Theorem 10), 



Q m {x;a,P, iV)) 



f) / ;/ \ n-m „| 

dl n [ ' ^{fM-l) m\(n-m) m [) 



^k%) (x, N) = (—1 + n — iV) (re, N) 



(see e.g. [12], Theorem 9) 

dp 

(see e.g. [12], Theorem 9), 



9 ft) ( X , N) = n (-1 + n - N) k^x, N) , 



a n—1 ji— m— 1 „i 

^ m=0 v ' 

o n—1 I 

—K^'flfx) = V - K^(x) 

8(3 n [X) ^ n a{n-m)m\ m [h 

Proof: If the derivative is taken with respect to a variable occurring as an argument rather 
than as a parameter in the hypergeometric representation, its representation can be easily 
obtained from the derivative rule of the generalized hypergeometric function, and the chain 
rule. In those cases, the representations need at most two neighboring polynomials. 
The other cases can be handled similarly to Corollary 3. □ 
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7 Conclusion 

Here, we want to recall the algorithms to convert between the different types of representations: 

1. Hypergeometric Representation — > Recurrence Equation: Zeilberger's algorithm 

2. Hypergeometric Representation — > Difference/Differential Equation: Zeil- 

berger's/Almkvist-Zeilberger's algorithm 

3. Difference/Differential Equation — > Recurrence Equation: Theorem 1 

4. Difference/Differential Equation — > Hypergeometric Representation: method 
of § 3 and § 5 

5. Recurrence Equation — > Difference/Differential Equation: Algorithms 1 and 2 
in [14] 

6. Recurrence Equation — > Hypergeometric Representation: combination of meth- 
ods 5 and 4 
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